Hopping conductivity of a suspension of nanowires in an insulator 
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We study the hopping conduction in a composite made of straight metallic nanowires randomly 
and isotropically suspended in an insulator. Uncontrolled donors and acceptors in the insulator lead 
to random charging of wires and, hence, to a finite bare density of states at the Fermi level. Then 
the Coulomb interactions between electrons of distant wires result in the soft Coulomb gap. At low 
temperatures the conductivity is due to variable range hopping of electrons between wires and obeys 
the Efros-Shklovskii (ES) law In cr oc -(T BS /T) 1/2 . We show that T E s oc l/(nL 3 ) 2 , where n is the 
concentration of wires and L is the wire length. Due to enhanced screening of Coulomb potentials, 
at large enough nL 3 the ES law is replaced by the Mott law. 
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In recent years, electron transport in composites made 
of metallic granules surrounded by some kind of insula- 
tor generated a lot of interest in both fundamental and 
applied research. Conductivities of the composites with 
spherical or generally speaking single scale granules are 
well studied (see review [if and references therein) . If the 
number of granules n per unit volume is large, granules 
touch each other and the conductivity is metallic. When 
n is smaller than the percolation threshold n c , granules 
are isolated from each other, the composite is on the in- 
sulating side of metal-insulator transition and the con- 
ductivity is due to hopping. It was observed that in this 
case, the temperature dependence of the conductivity a 
obeys 



<T = cr cM-(To/T) a }, 



(1) 



where To is a characteristic temperature and a = 1/2. 
This behavior of conductivity a was interpreted as the 
Efros-Shklovskii (ES) variable range hopping (VRH) be- 
tween granules [ij, |2| . The finite bare density of states 
at the Fermi level was attributed to the charged impuri- 
ties in the insulator playing the role of randomly biased 
gates. Then the interaction of electrons residing on dif- 
ferent dots creates the Coulomb gap, which leads to ES 
variable range hopping. 

For needle-like granules such as metallic whiskers or 
metallic nanotubes, the experimental situation is more 
complicated than for single scale ones. Recently ES tem- 
perature dependence of the conductivity was reported [3| 
for composites made of bundles of single- wall carbon nan- 
otubes (SWNT) suspended in insulating polymers. On 
the other hand, the Mott law with a = 1/4 was also ob- 
served in many nanotube based materials Q, |j, |a, |fj . In 
order to understand the puzzling crossover between ES 
and Mott hopping, in this paper, we study the low tem- 
perature ohmic dc transport in an isotropic suspension of 
elongated metallic granules in an insulating medium (see 
Fig. [Q. Below we call such a granule a wire for brevity. 

In this paper we restrict ourselves to scaling approxi- 
mation for the conductivity and to delineating the corre- 
sponding scaling regimes. In our scaling theory, we drop 




FIG. 1: A suspension of neutral metallic wires in an insulator 
with concentration 1/L 3 Cn<C 1/L 2 d. 



away both all numerical factors and, moreover, also all 
the logarithmic factors, which do exist in this problem, 
because it deals with strongly elongated cylinders. In 
this context, we will use symbol "=" to mean "equal up 
to a numerical coefficient or a logarithmic factor" . 

We assume each wire is a cylinder with the length L 
and the diameter d <C L. When the concentration of 
wires n is very small, the distance between them is much 
larger than L, we can easily apply the results obtained in 
the Ref. |2j . Here we consider a range of concentrations 
1/L 3 <C n <C 1/L 2 d, where n is so large that spheres 
built on each wire as the diameter strongly overlap, but 
the system is still far below the percolation threshold 
n c = 1/L 2 d. Percolation threshold at very small con- 
centration n c — 1/L 2 d is the result of large excluding 
volume L 2 d which one wire creates for centers of others. 
(At concentrations above 1/L 2 d it also becomes impos- 
sible to place wires randomly and isotropically because 
of nematic ordering [7| , but there is no such problem at 
n <C n c ). Our results are summarized in Fig. [3 in the 
plane of parameters T and nL 3 . The main result is that 
at low temperatures, we arrive at ES law 



a = ( 7 exp[-(r BS /T) 1/2 ] 



with 



T, 



ES 



(2) 



(3) 



Ka(nL 3 ) 2 ' 
where k is the dielectric constant and a is the tunneling 



decay length in the insulator. The large factor (nL 3 ) 2 
in Eq. is a result of the enhancement of both the 
effective dielectric constant and the localization length 
in the composite. At large enough nL 3 and higher tem- 
peratures, ES law is replaced by the Mott law. And 
when the temperatures are sufficiently high, both VRH 
regimes are replaced by the activated nearest neighboring 
hopping (NNH) regimes. 

Let us start from the dielectric constant of the compos- 
ite. In the range of concentrations 1/L 3 <C n <C 1/L 2 d, 
its macroscopic dielectric constant n e s is greatly en- 
hanced due to polarization of long metallic wires. It was 
shown in the Ref 8, 91, that 



K e B = {nL 3 )n. 



(4) 



Such result can be understood in the following way. If 
the wave vector q is larger than 1/L, the static dielectric 
function for the wire suspension has the metal-like form 



e(q) = K 1 



1 



2 2 

q ri 



(5) 



where r s = (nL)~ l > 2 is the typical separation of the 
wire from other wires. It plays the role of screening ra- 
dius for the wire charge. The function e(q) grows with 
decreasing q until q = 1/L where the composite loses 
its metallic response and e(q) saturates. As a result, 
the macroscopic effective dielectric constant is given by 
K e fi = ({q = 1/L) = nL 3 K. Within our scaling approx- 
imation, one can also derive Eq. J2J starting from the 
facts that each isolated wire has polarizability equal to 
L 3 and due to random positions of wires the acting elec- 
tric field does not differ from the average one. 

The spacing between charging energy levels of such 
wire is the order of E c = e 2 /k c hL = e 2 / ' nnL A . If the wire 
is thick enough, it can be treated as a three dimensional 
object and the mean quantum level spacing 5 in the wire 
is the order of Xph /mLd 2 , where Xf is the Fermi wave- 
length and m is the effective electron mass. Narrow wires 
such as SWNT with single conducting channel will be dis- 
cussed later. Using the tunneling decay length a in the 
insulator, we can rewrite S in the form S = f3e 2 a 2 /nd 2 L 
where /3 = Xp/a. The ratio E c /S decreases proportional 
to 1/nL 3 because of enhanced screening of Coulomb in- 
teraction in the denser system of longer wires. 

We start our discussion with relatively small nL 3 such 
that nL 3 < d 2 j (5 a 2 and E c > 8. The density of states 
and localization length are all the information we need to 
calculate the VRH conductivity. In contrary to a doped 
semiconductor, in a suspension of neutral wires, the bare 
density of states at the Fermi level go — 0. Charging 
energy levels of the suspended wires are shown in Fig. 
02 a), they are equally spaced by E c . The Fermi level is 
at zero energy between empty and filled levels. 

Since the ground states of wires determine the low tem- 
perature hopping transport, we consider only the ground 
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FIG. 2: (Color online) The summary of the transport regimes 
for relatively clean wires schematically plotted in the plane 
of nL 3 and temperature T. The dark line separates regimes 
where Coulomb interaction is important (Ni [Eq. @] and ES 
law [Eq. J3J]) from regimes where it plays a minor role (Na 
[ ifTHjl ] and Mott (M) law [Eq. E|). The grey (blue) line sepa- 
rates activated the nearest neighbor hopping (NNH) regimes 
from the VRH regimes. Shaded domain represents the metal- 
lic side of the insulator-metal transition. Instead of smooth 
crossovers at dark and grey lines, in the vicinity of this tran- 
sition, the conductivity has a critical behavior. 
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FIG. 3: The density of ground states at E c > S. (a) BDOGS 
in neutral wires, (b) BDOGS for charged wires with large 
enough concentration of donors in the insulator in the absence 
of the long range Coulomb interaction; (c) The density of 
ground states of charged wires with the Coulomb gap in the 
vicinity of the Fermi level fi. Occupied states are shaded. 



state in each wire at a given number of electrons and 
exclude excited states 2J. Thus, the density of states 
we need can be called bare density of ground states 
(BDOGS). According to Ref. 0, the finite BDOGS near 
the Fermi level originates from uncontrollable donors (or 
acceptors) in the insulating host. Donors have the elec- 
tron energy above the Fermi energy of wires. Therefore, 
they donate electrons to wires. A positively charged ion- 
ized donor can attract and effectively bind fractional neg- 
ative charges on all neighboring wires, leaving the rest of 
each wire fractionally charged. At a large enough aver- 
age number (3> 1) of donors per wire, effective fractional 
charges on different wires are uniformly distributed from 
— e/2 to e/2. In such a way the Coulomb blockade in a 
single wire is lifted and the discrete BDOGS get smeared 
(see Fig. El»). The BDOGS g Q becomes 1/L C L 3 . In 
the very vicinity of the Fermi energy, the long range 
Coulomb interaction creates the parabolic Coulomb gap 



A (see Fig. 0{ c ))- The Coulomb gap affects energy in- 
terval |e| < A. At low enough T, when the range of 
energy levels around the Fermi level responsible for hop- 
ping {TesT) 1 / 2 is much smaller than A, the conductivity 
ob eys the ES law (see Eq. J3J)) and does not depend on 
go [l0(- The parameter Tes — e 2 /K e s£,, where £ is the 
localization length for tunneling to distances much larger 
than L. Let us concentrate now on the nontrivial value 
of the localization length £. 

It was suggested [llj that such long range hopping pro- 
cess can be realized by tunneling through a sequence 
of wires (see Fig. 0J. The states of the intermedi- 
ate wires participate in the tunneling process as virtual 
states. Nowadays the virtual electron tunneling through 
a single granule is called co-tunneling [lj, |l2, LLiJ LL-I \lli 
and regarded as a key mechanism of low temperature 
charge transport via quantum dots. One should distin- 
guish the two co-tunneling mechanisms, elastic and in- 
elastic. During the process of elastic co-tunncling, the 
electron tunneling through an intermediate virtual state 
in the granule leaves the granule with the same energy as 
its initial state. On the contrary, the tunneling electron 
in the inelastic co-tunneling mechanism leaves the gran- 
ule with an excited electron- hole pair behind it. Which 
mechanism dominates the transport depends on the tem- 
perature. Inelastic co-tunneling dominates at T > \/E c 5, 
while below this temperature, elastic co-tunneling wins. 
In discussions below, we consider sufficiently low temper- 
atures such that only elastic co-tunneling is involved. 

First we remind the idea of calculation of the localiza- 
tion length £ used in R.e f. 11 ll . (Such approach was also 
applied in the papers |2,llj|). When an electron tunnels 
through the insulator between wires at the nearest neigh- 
boring hopping (NNH) distance r, it accumulates dimen- 
sionless action r/a, where a is the tunneling decay length 
in the insulator. We want to emphasize that the NNH 
distance r is realized only in one point of the wire and, 
therefore, is much shorter than the typical separation 
along the wire from other wires r s . NNH distance r can 
be calculated using the percolation method |lO| . If one 
builds around each wire a cylinder with length L and ra- 
dius r, percolation through these cylinders appears, when 
nL 2 r — 1. This happens because excluded volume cre- 
ated by one cylinder is not r 2 L, but L 2 r. As a result, 
we obtain r — 1/nL 2 . We assume the metallic wire is 
only weakly disordered and hence we neglect the decay 
of electron wave function in the wire. (This assumption 
will be relaxed later.) Over distance x, electron accu- 
mulates x/L additive actions of the order of r/a. Thus, 
its wave function decays as exp(— xr/La) = exp(— x/£), 
where the localization length £ = aL/r = a(nL 3 ). This 
contributes another factor nL 3 to Eq. @. Eq. (0) with 
Tes given by Eq. © is valid in the left lower corner of 
Fig. ^labeled as ES. More rigorous calculation [l£j of the 
localization length £ gives £ = L/[ln(E c ir/5) + l/nL 2 a\. 
Our simple derivation of £ along the line of Ref. [ll| has 




FIG. 4: (Color online) An illustration of a long range (2> L) 
hop from the wire 1 to the wire 2 via a sequence of wires. 



lost the first term in the denominator. We argue that in 
our case l/nL 2 a ^> d/a, and d/a can be quite large. As 
a result, within a large range of the ratio E c /S, the term 
ln(E c Tr5) is only a small correction to the leading term 
l/nL 2 a. 

At high temperatures, the conductivity is dominated 
by the nearest neighbor hopping. Using the NNH dis- 
tance r, we obtain the conductivity with activated T- 
dependence, 



<Toce 



-1/nL 2 



-E A /T 



(0) 



where the activation energy Ea is determined by the 
charging energy 



E A = 



k cS L 



iL* 



if 5-€.E c 



(7) 



The range of validity of Eqs. JJjJ and JJJ) is shown in the 
upper left corner of Fig. |21 and labeled as Ni . 

When L becomes so large that S > E c , the BDOGS go 
evolves from {E C L 3 )- 1 to g = {SL 3 )- 1 = nd 2 /e 2 pL 2 a 2 . 
Since go decreases and k c s grows with L, the width of the 
Coulomb gap A, which can be estimated as = y^ goe e / n 3 s 
decreases with L. Eventually, the Coulomb gap becomes 
narrower than the width of Mott's optimal band €m = 
T 3/4 /(gof) 1/4: at given T. At this point, the ES law is 
replaced by the conventional Mott's law, 



a exp[-(r M /r) 1 / 4 ], 



where 



T M = 



(3e 2 



got 3 



i 3 L 7 d 2 ' 



(») 



(9) 



This regime is labeled as M in Fig. [5] It crossover to the 
ES regime at T = e 2 d 2 / Ka(3nL 5 , which can be obtained 
by equating A and em- This change of the temperature 
dependence of the conductivity from 1/2 ES law to 1/4 
Mott law was observed at relatively high temperatures 
and large concentrations in the experiment [3|- 

At high temperatures, the Mott's law is replaced by 
the NNH conductivity with activation energy equal to 
the spacing of quantum levels 8. The T-dependence of 
the conductivity is, 
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FIG. 5: (Color online) Transport regimes for wires, where 
electrons are localized at length X < L. Dark line separates 
regimes where Coulomb interaction is important (Ni_3 and ES 
laws) from regimes where it plays minor role (N2,4 and Mott 
laws). The grey (blue) line separates activated NNH regimes 
from VRH regimes. The left dashed line marks the border 
where dielectric constant K c g starts growing with nL 3 , while 
the right dashed line serves as the border where localization 
length £ stops growing and becomes A. The shaded domain 
represents the metallic side of the insulator-metal transition. 
Instead of smooth crossover, in the vicinity of the transition, 
the conductivity has a critical behavior. 



We call this regime N2 and show it in the right upper cor- 
ner of Fig. It is easy to check that ES and Mott VRH 
transport regimes match corresponding NNH regimes Ni 
and N2 smoothly, when VRH distance is of the order of 
the wire length L. This happens at the gray line of Fig. 
|2| The Coulomb regimes ES and Ni smoothly crossover 
to the regimes M and N2 of noninteracting electrons at 
the dark line of Fig. [21 We stop at nL 3 = L/d, where 
percolation leads to the insulator-metal transition. Our 
scaling approach is not designed to study the critical be- 
havior of the conductivity in the vicinity of the transition. 
Now we switch to the case of metallic wires with a 
strong internal disorder, where the electron wave func- 
tion decay length A is much smaller than the wire length 
L. The summary of transport regimes for this case is 
presented in Fig. [S] In order to show all the possible 
regimes, we assume that A is small enough for the order 
of labels on the nL 3 axis shown in Fig. |SJ (At a larger A 
only some of the shown regimes may exist.) We also plot- 
ted regimes labeled as ES and Ni on the left of the border 



line nL 



3 _ 



iX . The characteristic energies of electron 



transport in these regimes are given by Eqs. 10 and Q 
respectively. When L ^S> A, the wires can be treated as a 
larger number of clean metallic wires with length of the 
order A and the concentration nL/X. Following the logic 
of Ref. [8j and our derivation of n g (Eq. (J3J)) one can 
show that the effective macroscopic dielectric constant is 
given by 



^offl 



kjiLA . 



(11) 



Comparing to (Eq. Q) we see that the dielectric con- 
stant K e ffi is strongly reduced by the internal disorder 
of wires. This happens because the disorder effectively 



"cuts" the wire into short pieces albeit increases the con- 
centration of the short wires. 

The internal disorder also changes the localization 
length for the long distance (3> L) tunneling. Using the 
argument we made above, we obtain 



1 



iL 3 a 1 



(12) 



where the first term on the right hand side accounts for 
the decay of the wave function in the wire, while the sec- 
ond term represents the decay in the insulating medium 
between the wires. When nL 3 < X/a, the second term 
dominates, £ = nL 3 a, which is the same as what we got 
for the case of clean wire. However, because of new di- 
electric constant K e si (Eq. Illfl . we obtain new regime 
labeled ESi in the Fig. \5\ Here the conductivity obeys 
ES law with the characteristic temperature 



Tes x — 



Kan 2 L 4 X 2 



(13) 



The charging energy E c also changes to e 2 / 'K e siL = 
e 2 / nnL 2 X 2 . It decreases with growing nL 3 and reaches 
the mean quantum level spacing 5 of the wire at nL 3 — 
d 2 L 2 / j3a 2 X 2 . With relatively high temperature and small 
nL 3 < d 2 L 2 /j3a 2 X 2 , the regime ESi crossover to the 
regime N 3 , where the conductivity can be represented 
by the same equation © but with a different activation 
energy Ea — e 2 /nnL 2 X 2 . When nL 3 > d 2 L 2 / (3a 2 X 2 , the 
Coulomb interaction plays a minor role. Hopping con- 
ductivities in regimes M and N2 are not affected by the 
change of the dielectric constant and are given by the 
same Eqs. JHJ, © and ©, (fTTTft respectively. 

At even larger nL 3 > X/a, the decay of electron wave 
functions in the wire dominates, thus £ = A. As a result, 
we have another ES regime, ES2, with 



Tes 2 = 



the Mott regime Mi with 



1LX 3 ' 



Tm 1 = 



i 2 L 4 X 2 ' 



(14) 



(15) 



and the activated nearest neighboring hopping regime 
N4, with the activation energy S and 



a oc exp l — — ] exp 



e 2 Pa 2 \ 
' Kd 2 LTJ 



(16) 



Thus, we have completed our consideration of the 
phase diagram up to nL 3 = L/d, at which wires begin to 
percolate. What happens to the conductivity if one man- 
ages to create an isotropic suspension with nL 3 > L/d? 
In contrary to relatively clean wires, there is no insulator- 
metal transition around nL 3 = L/d, where wires start to 
touch each other. Indeed such transition happens when 



the typical branch of the percolating network is metallic, 
which means the typical distance between two nearest 
neighboring contacts along the same wire can not be L, 
but should be smaller than A. We can now think about 
percolation over pieces of wires with the length A and 
the effective concentration nL/X. Such percolation ap- 
pears at (nL/\)(\ 2 d) = 1 or nL 3 = L 2 /\d. This is the 
threshold of insulator-metal transition. As a result our 
phase diagram of hopping regimes is extended as far as 
nL 3 = L 2 /Xd. 

In this paper we concentrated on relatively thick wires 
with the three dimensional density of states. Narrow 
wires with as few as one conducting channel should 
be treated as one dimensional objects. In such case 
the level spacing is given by 5 — fi /ttiXfL and thus 
E c /S = (l/nL 3 )(X F /a) — /3/nL 3 . Since most probably 
(3 is the order of 1, while nL 3 ^> 1, the charging energy 
E c is always smaller than level spacing S. Coulomb in- 
teraction plays minor role in such system. Therefore the 
left part of our phase diagram Fig. ^including ES regime 
vanishes. 

It is tempting to apply this theory to carbon nan- 
otubes. Metallic SWNT has only one conducting chan- 
nel. Therefore for them, the domain of ES hoping most 
likely vanishes. Bundles of SWNT, studied in Ref. |g def- 
initely have many channels and our theory should work 
for them providing the geometry of bundles may be ap- 
proximated by a set of randomly oriented rigid rods of 
a given characteristic length L. In this case, the range 
of ES conductivity should exist. ES range may also ex- 
ist for multi-wall nanotubes, where one can expect many 
parallel metallic channels. 

In this paper, we considered only straight wires. One 
can generalize this theory to random system of flexi- 
ble wires formed by interpenetrating Gaussian coils (for 



example conducting polymers suspended in insulating 
medium). This theory will be published elsewhere. 
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